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Within the Dyson-Maleev boson formahsm, we study the zero-temperature magnon dispersions 
in a family of one-dimensional quantum Heisenberg ferrimagnets composed of two different spins 
(Si ,5*2) in the elementary cell. It is shown that the spin-wave theory can produce precise quantitative 
results for the low-energy excitations. The spin-stiffness constant ps and the optical magnon gap A of 
different (S\, S2) ferrimagnetic systems are calculated, respectively, to second and third order in the 
quasiparticle interaction. The spin-wave results are compared with available numerical estimates. 
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I. INTRODUCTION 

In the last decade a large variety of quasi-one- 
dimensional (ID) mixed-spin compounds sMith ferri- 
magnetic properties has been synthesized. B Most of 
them are molecular magnets containing two different 
transition-metal magnetic ions which are alternatively 
distributed on the lattice (Fig. |l|). For example, two 
families of such compounds read ACu(pba)(H20)3-nH20 
and ACu(pbaOH)(H20)3-nH20, where pba = 1,3- 
propylenebis (oxamato), pbaOH = 2-hydroxo-l,3- 
propylenebis (oxamato), and A=Ni, Fe, Co, and Mn. 
Published experimental work implies that the magnetic 
properties of these mixed-spin materials are basically 
described by a quantum Heisenberg spin model with 
antiferromagnetically coupled nearest-neighbor localized 
spins: 

N 

n - [Si(n) + Si(n + 1)] • S2(n) 

n=l 

N 

- lisH J2 [giS^H + g2SI{n)] , J > . (1) 



Here the integers n number the N elementary cells, each 
of them containing two site spins with quantum num- 
bers S'l > 82- The 5-factors, related to the spins 
and S2, are denoted by gi and 32, respectively. fiB 
is the Bohr magneton, and H is the external uniform 
magnetic field applied along the z direction. As an ex- 
ample, the following values for the parameters in Eq. 
(|l|) have been extracted from magnetic measurements on 
the recently synthesized quasi-lD bimetallic compound 
NiCu(pba)(D20)3-2D20: (51,^2) = {Sm,Scu) - (Ipi), 
J/kB = 121K, g, = gN, = 2.22, 32 = gcu = 2.09.1 In 
view of the recent developments in the physics of uniform 
quasi-lD systems (concerning, in particular, the behav- 
ior of antiferromagnetic chains and ladders in external 
magnetic fields: see, e.g., Ref. ^ and references therein), 
ferrimagnetic chains and ladders open an interesting new 



area. The point is that the presence of two or more dif- 
ferent quantum spins in the elementary cell considerably 
increases the number of situations of interest: the topol- 
ogy of spin arrangements plays an essential role in tha 
structure of the ground state and low-energy excitations.cl 
For example, the diamond lattice in Fig. |l| represents an- 
other interesting class of Heisenberg ferrimagnetic chains 
constructed of one kind of antiferromagnetically coupled 
site spins but with a dificrent number of lattice sites in 
the A and B sublatticesE It is remarkable that only in 
the last few years interest in the physical properties of 
ID quantum Heisenberg ferrimagnets has considerably 
increased:^ most of the early efforts were concentrated on 
the chemistry of molecular magnets and a relatively sxnall 
amount of work was devoted to physical properties .0 
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FIG. 1. Two types of ferrimagnetic chains on bipartite 
lattices: a) mixed-spin chain with alternating site spins 5*1 
and S2 {Si > S2)', b) uniform-spin diamond chain with a site 
spin S. 

Since interactions in the models presented in Fig. |^ 
set up bipartite lattices, the Lieb-Mattis theoreml is 
applicable. In particular, for the mixed-spin system it 
predicts the existence of a ferrimagnetic ground state 
with the total-spin quantum number Sg = {Si — S2)N. 
Thus the model (|l|) exhibits long-range ordered magnetic 
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ground state characterized by the magnetization density 
Mq = (^i — 52)7200, ao being the lattice constant. Such 
a magnetic phase may be referred to as a quantized unsat- 
urated ferromagnetic phase: it is characterized by both 
the quantized ferromagnetic order parameter M (quan- 
tized in integral or half-integral multiples of the num- 
ber of elementary cells, M = MqN) and the macro- 
scopic sublattice magnetizations Ma = X^„=i(Si(?t.)) 

and M.B ~ J2n=ii^^(''^)) ■ What makes such a state and 
the related transitions to magnetically disordered states 
interesting is the fact that the ferromagnetic order pa- 
rameter is a conserved quantity: in Heisenberg systems 
this is expected to imply strong constraints on the criti- 
cal behavior.u Consequences of the spontaneous symme- 
try breaking SO{3) — > 50(2) (related to the establish- 
ment of a ferrimagnetic ground state) for the structure of 
low-energy excitations are dictated by the nonrelativistic 
version of Goldstone's theoremll3: in the absence of long- 
range forces, a spontaneous symmetry breaking leads to 
low-energy excitations whose energy tends to zero for 
wave-vectors fc — > 0. In contrast to the relativistic ver- 
sion, the theorem does neither specify the exact form of 
the dispersion relation for small fc, nor does it determine 
the number of different Goldstone modes: these features 
are not fixed by symmetry considerations alone ~ rather, 
they depend on the specific properties of the system. 

Turning to the case of Heisenberg ferromagnets, in 
a hydrodynamic description the quadratic form of the 
magnon energies 

E,^-^^k^ + 0{k^) (2) 

results entirely from the symmetry of the ferromagnetic 
state and the fact thatj-the order parameter is itself a 
constant of the motion.Ell Here Mq and ps are, respec- 
tively, the magnetization density and the spin-stiffness 
constant of the Heisenberg ferromagnet. This form of 
the Goldstone modes (which is just the Landau-Lifshitz 
formula) can rigorously be obtained by simple sum-rule 
arguments (see, e.g., Ref. |l^). Similar argiimpnts are 
applicable to Heisenberg ferrimagnets as well.EjLJ An al- 
ternative approach, relying on the suggested conformal 
invariance of the^odel also predicts the quadratic 
form of Eq. (||)Jl3 Due to the doubling of elementary 
cell in the ferrimagnetic ground state, there appears a 
second spin-wave branch in the ferrimagnetic case (op- 
tical magnons) with a finite energy gap A at fc = 0: 
A ^ Afo/x||,|-iyhere x\\ is the magnetic susceptibility 
parallel to M.ll^ At intermediate temperatures, it is the 
optical magnon branch which produces a number of spe- 
cific thermodynamic properties-characteristic especially 
for the ferrimagnetic systems. □'EZTlB 

The existence of a macroscopic magnetic ground 
state also opens an interesting and rare opportunity 
to apply the spin-wave theory (SWT) to the low- 
dimensional quantum spin system (|l|). Indeed, re- 
cently published ai aa J ^r s is . based on the linear spin- 
wave approximationOEllES qualitatively confirmed the 



expected structure of low-energy excitations in Heisen- 
berg ferrirpagnets. Moreover, the second-order SWT 
was shownE^ to produce precise quantitative results for 
the ground-state energy Eq and the on-site magnetiza- 
tions nil — Ma/N and m2 = Mb/N even for the ex- 
treme quantum system {81,82) — (1,5): the second- 
order SWT results differ by less than 0.017% for Eq, and 
by less than 0.18% for mi , from the density matrix renor- 
malization group (DMRG) estimatesJlZl In this respect, 
an interesting question is to what extend the spin-wave 
approach can effectively be used to describe the proper- 
ties of the ID model (|l|) at zero temperature. A purpose 
of the present paper is to demonstrate, through an ex- 
plicit study of the magnon-dispersion perturbation series, 
that the spin-wave approach can produce precise quanti- 
tative results for the magnon dispersions as well. 

The paper is organized as follows. In Section 2, using 
the Dayson-Maleev boson representation of spin oper- 
ators, the original spin Hamiltonian (|^) is transformed 
to an equivalent boson Hamiltonian. The choice of an 
appropriate zeroth-order quadratic Hamiltonian for the 
perturbation series is also discussed here. In Section 3 
we study the second-order corrections for the magnon 
energies and the related corrections for the spin-stiffness 
constant and the optical magnon gap. Third-order self- 
energy diagrams and the respective corrections for the 
optical magnon gap are also considered here. In the 
last Section the results are summarized and a compar- 
ison with available numerical estimates is made. 



II. DYSON-MALEEV FORMALISM 

To develop a spin- wave theory (see, e.g., Refs. |24|j25| 
and references therein), one firstly transforms the origi- 
nal spin Hamiltonian (n|) to an equivalent boson Hamil- 
tonian. We adopt the following Dyson-Maleev represen- 
tation of spin operators: 

8^ {i) = \/2S^ (^ai - '^'^i^i^i^ ' (3) 

5f (z) - y^ial , 8f {t) = 8i-ala,, (4) 

8+{j) = (^b] - ^/A^j^ , (5) 

8^{j)^V^2b„ 8^{j)^b%~82, (6) 

where 8^ = 8^ ± i8]^, a = 1,2. ai and bj are boson 
operators defined on the lattice sites i G A and j E B, 
respectively. 

Substituting the latter expressions in Eq. (^ , one finds 
the following boson representation of TL in terms of the 
Fourier transforms and b^ of a.i and bj : 

He = -2c78^NJ + n'o + V^m , (7) 

where 



2 



and 



V 



DM 



N 



(8) 



(9) 



( 2-fi^4ala2b\b4 + y/aji+2-4:alblblbi + -^7403020164). 

7fc — cos(fcoo) is the lattice structure factor and = 
(5(1 + 2 — 3 — 4) is the Kronecker 5- function. We have 
used the convention (^1,^2,^3,^4) = (1,2,3,4) and the 
notations a = S1/S2 > 1 and 5*2 = S. Here and in what 
foUows the external magnetic field is 7J = 0. 

In the above expressions Hq — 0{S) is the quadratic 
boson Hamiltonian of the linear spin-wave theory 
(LSWT) and V^]^.i ^ 0{1) is the quartic Dyson-Maleev 
boson interaction. Omitting completely the boson in- 
teraction Vjjj^j, a number of authors has recently used 
the quadratic LSWT Hamiltonian TYq for a qualitative 
description o£n th| i i) .ID ferrimagnetic model (yj) at zero 
temperatureJl3'EilE3 Tig is easily diagonalized by use of 
the following Bogoliubov transformation to the quasipar- 
ticle boson operators ak and Pk'- 



Ofe = Ukiak 



Uk 



l + Sk 



2ek 



Xk 



bk = UkiPk 
1 + efe ' 



XkCtl) , 



(10) 



where Sk = \/l-r]k ^nd r]k = 2^/^7^/(1 + a). 

As a matter of fact, for the ID ferrimagnetic model (|^) 
the quadratic Hamiltonian (^) is not the most appropri- 
ate choice of a starting zeroth-order approximation for 
perturbation series. Since we are also interested in sys- 
tems with small site spins, it is more appropriate from 
the very beginning to redefine Hq by addiM the first- 
order 1/5 corrections to magnon dispersions .Ej The latter 
corrections come entirely from a normal ordering of the 
quasiparticle boson operators ak and (3k in the quartic 
interaction V^j^j. It is important that these corrections 
(similar to Oguchi's corrections in antiferromagnetsEZI) 
renormalize the magnon excitation spectra (and the 
ground-state energy) without changing their basic struc- 
ture, i.e., the number of Goldstone modes. Unlike the 
case of Heisenberg antiferromagncts, where Oguchi's cor- 
rections are numerically small even for site spins S* = ^, 
in quantum ID ferrimagnets the magnon dispersions are 
significantly renormalizedcj (see below). As a result of 
the normal-ordering procedure, the boson Hamiltonian 
(^ is recast to the following basic form: 



Eo+Hn + XV, V = V2 + Vi 



DM 



A = 1, 



(11) 



The ground-state energy Eq calculated up to first order 
in l/S* reads 



Eo 
2NJ 



= -aS' 



S [2y^ai + (cr + 1)02] 
cr + 1 



-ai02 , 



(12) 



where 



Oi 



2(a+l)7V^£fc' ' 2 2N^ek' 

(13) 



Expressed in terms of quasiparticle operators, the cor- 
rected LSWT Hamiltonian reads 



no^2SjJ2[uji:'^alak+4^^(3i(3k 



(14) 



where 



01 \ /tr + l a-l\ 1-72 
: I I — ■ — ei- T — ■ — I - 02 ^ ^ 



(15) 



The normal-ordered quasiparticle interaction V con- 
tains a quadratic term. 



where 



V2 



(+,-) 



1± 



is.) 2^ 

1 J Sk 



(16) 



(17) 



and the quartic normal-ordered Dyson-Maleev interac- 
tion Vdm containing nine vertex functions: F*^*-* = T^i'2^34, 
i — 1, ... ,9. Exphcit expressions for Vdm and the related 
vertex functions are presented in Appendix A. 

The quadratic Hamiltonian (|l^) describes two 
branches of spin- wave excitations (magnons) defined with 



the dispersion relations E 



(",/3) 



= 2SJuj 



Eq.(|l|). 

The a excitations are gapless acoustical magnons in the 
subspace with Sf^^ = {Si — S2)N — 1, whereas the /3 
excitations are gapful optical magnons in the subspace 
S^,, = (5i - S2)N + 1: Stot = ELi [Si(") + S2(n)]. It 
is easy to show that the effect of the external magnetic 
field in Eq. (|^) on the dispersions eI,"'^^ reduces to the 
simple substitution e\^'^'^ \ — > e'^'^'' ±g^BH provided 
that 51 = ^2 = 9- 

The spin stiffness constant ps (related to the acoustical 
branch) plays a basic role (together with the magnetiza- 
tion density Mq) in |y:u3, lejaihtemperature thermodynam- 
ics of the model (|^)la't3'E§Ej and can be obtained from 
the Landau-Lifshitz relation (||) and Eq. (|l5|): 



J Oq S'iS'2 



ai 



1 



02- 



Sa 



(18) 
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In the same free-quasiparticle approximaticp, the spec- 
tral gap of optical magnons at k — readfES 



Ao = 2J(5i -^2)1 



(19) 



and (|T|) give 
On the other 



In the case {81,82) = (1, 5) Eqs. (|18|) 
p's^^/Jao8i82 = 0.761 and Aq = 1.676J 
hand, the LSWT Hamiltonian Hq produces the parame- 
ters of the related classical system: pi°^/JaoS'iS'2 = 1 
and Aq = pi. The numerical estimate for the gap 
A = 1.759jE2l clearly demonstrates the importance of 
the 1/8 corrections in Eq. ([l5|). 

Next, let us consider the macroscopic sublattice 
magnetizationsQ mi = {8l{i)) — 81 — {o\a,i) and m2 = 
(S'l(j)) = —5*2 -I- {b^jbj) which are finite in the ferrimag- 
netic state. In the case of small site spins and for small 
magnetization densities Mg, mi and m2 are strongly re- 
duced as compared to their classical values (respectively, 
and —82)'- 



nil = 81 - a2 , mi + m2 = 5i — S'2 . 



(20) 



As an example, for (^i, 82) — (1, ^) Eq. (20) predicts a 
42% reduction of the small spin 82 — ^- However, the 
off-diagonal quadratic interaction V2 in Eq. (|l|) produces 
important first-order corrections for the sublattice mag- 
netizations mi and m2. Thus, to first order in 1/8 Eq. 
( pO| ) should be replaced by the more precise expression 



mi = ^i - a2 - 



1 



28{a + l)^N^ 



r-\ Ik 



vn^- (21) 



Now for the ferrimagnetic chain (1,5) the last expression 
gives the result mi = 0.816 which diffcrs_by only 3% from 
the DMRG numerical estimate 0.79248tZl. 

Summarizing, it may be stated that the free- 
quasiparticle approximation based on the Hamiltonian 
(1^) gives a good qualitative description of the ground- 
state properties of the model (|l|). Further improvement 
of the SWT results may be achieved by considering the 
role of quasiparticle interactions. 



III. ROLE OF THE QUASIPARTICLE 
INTERACTIONS 

Since the first-order 1/8 corrections for the ground 

state energy Eq and the dispersions E^''^^ have already 
been taken into account by the normal-ordering proce- 
dure, corrections to Eqs. (|l^) and (|l^) arise only up 
from the second order of the perturbation series in V. 



the self-energy diagrams in Fig. 
lytic expressions read 



The respective ana- 



<54"'(a) = 



(+)t/(-) 



vr'v, 



{28 f J^) 



(22) 



J- ^ ^kp-pk ^ ^kp-pk 

2^ ..(") , ,.(/5) ' ^^'^> 



{28 f N 



{28f 



^34 



^(8) y(l) 
•^43:2*: *^fc2;34 



(24) 



-5-1"' (e) = 



{2Sf iV2 



34 

fc2" 



y(3) T.(2) 
'^43;2fe ''fe2;34 



-LO 



(a) 



(25) 



(3) 



Note that since the vertex functions Vj:^^ Vkp-pk^ ^kp-pk^ 
vj^^k^ ^^'^ ^4*3-2*: vanish at the zone center k = (see 
Appendix A), the gapless structure of w^"' is preserved 
by each of the second-order corrections, Eqs. (p^) to ([2^). 



v(+) v(-) 

< » <- 




a) 



< • > f< 




A. Second-order corrections for the magnon energies 

The second-order corrections for the dispersion of 
acoustical magnons w^"', Eq. ([l5|), are connected with 



FIG. 2. Second-order self-energy diagrams giving the cor- 
rections for the dispersion of acoustical magnons Solid 
and dashed lines represent, respectively, the propagators of 
a and /3 magnons. The Dyson-Maleev vertex functions V^^\ 
i = 1, . . . ,9 are defined in Appendix A. 
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The second-order corrections for the dispersion of opti- 
cal magnons io\f\ Eq. (|l5|), come from similar diagrams 

(see Fig. ^). The first diagram gives Sujlf \a) = (5a;["\a), 
whereas the other four diagrams give the following con- 
tributions: 



1 



{2Sy N 



E 



kp;pk 



(-)t/(6) 



P 



kp;pk 



(26) 



y(7} yiS) 



(25)2 ^2 /c2^(^) 



-|- UJ. 



(a) 



(2^)2 iV2 



,(/3) 



v(+) 



v(+) v<-) 
■-■>• > •■ > 



(27) 



^(5) y(6) 
•^43:2*; ^fc2;34 



,(/3) 




(28) 




Oguchi's treatment of the Heisenberg antiferromagne 
and is supported by the following observations: (i) the 
higher-order corrections to the principle approximation 



for 



,(a,/3) 



are numerically small (see below), and (ii) the 



third-order series in V gives a somewhat better result for 
the gap A in the extreme quantum system (1, ^)- As a 
matter of fact, noticeable deviations from the standard 
1/5 expansions appear only in the above mentioned ex- 
treme quantum cases. 

The dispersion functions E^'^^ for the system (1, ^), 
calcutated up to second order in V, are presented in 
Fig. IJ. For the optical magnon branch we find an excek 
lent agreement with the exact-diagonalization resultsE2l 
in the whole Brillouin zone. Recently, a very successful 
description of the optical magnon branch of the system 
(1, i) has also been achieved through the variational ma- 
trix product approach.E3. As to the acoustical branch, 
the agreement with the available quantum Monte Carlo 
resulta23 is not so satisfactory, especially near the zone 
boundary at fc = 7r/2ao. It is worth noticing that in the 
above region the calculated second-order corrections to 
ui^\ Eq. ([l5|), are very small (about 0.5% of the prin- 
cipal approximation). In principle, it is not excluded 
that the spin-wave expansion partially fails to describe 
the acoustical branch in the system (1, ^), as the latter 
is characterized by the minimal magnetization density 
Mq — l/4ao (see below). To clearify the problem, fur- 
ther studies are needed in this direction. 
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FIG. 3. Second-order self-energy diagrams giving the cor- 



rections for the dispersion of optical magnons ui^f^ (see also 
the notations in Fig. ^). 
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Note that in the present perturbation scheme (u sing 
Tio as a zeroth-order Hamiltonian) , in Eqs. (E^) to ( |2^ ) 
there appear the renormalized dispersions lo^'^\ The 
standard 0{1/ S"^) corrections to the magnon dispersions 
can easily be obtained by substituting in Eqs. ( p2[ ) to 

( p8| ) the bare excitation energies (i.e., the functions u^"''^'' 
without the 1/S corrections). Since we are also inter- 
ested in the extreme quantum systems which are com- 
posed of small site spins and have small magnetization 
densities Mq [such as (1, 5) and (|,1)], it may be ex- 
pected that the adopted perturbation scheme, where the 
quasiparticle interaction V is treated as a small perturba- 
tion, is more appropriate. Such a viewpoint is similar to 



FIG 



4. Dispersions of the acoustical S^"' and optical E^J^^ 



magnons (in J units) in the system {81^82) = (1, |) calcu- 
lated up to second order in the magnon interaction V. The 
points (-1- and o) denote, respectively, quantum Monte Carlo 
and exact-diagonalization results (Ref. |30[ ). 

The above second-order corrections to ^^^'^^ can be 
used to find the coefficient r2 in the perturbation series 
for the spin-stiffness constant 



J ao S1S2 



2S 



r2 



{2SY 



o 



(29) 
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and the coefficient S2 in the series for the optical magnon 
gap 



B. Third-order corrections for the optical magnon 

gap 



A 



2J(5i - ^2) 



Si S2 ^3 
2S ^ J2S)^ ^ J2S)^ 



O 



(30) 



Here ri — ~1a\l yfa — la^ia + 1)/(t and S\ — ~1a\i 
are obtained from Eqs. ( pq ) and (p9|), respectively. The 
results for r2 and ^2 are presented in Table 1. 



The third-order corrections for the optical magnon 
gap A are connected with the self-energy diagrams in 
Fig. § These can be obtained by drawing all connected 
irreducible self-energy diagrams containing three vertex 
functions and only oppositely-oriented (a, /3) pairs of in- 
ternal magnon lines. The diagrams containing magnon 
lines closing on themselves vanish since the vertices have 
already been normal ordered. In addition, the diagrams 
which have internal a lines carrying the momentum fc 
give zero contributions as well, since all vertex functions 
containing in-going a lines with momentum k vanish at 
fc = (see Appendix A). 




65) eg) 67) fig) 



FIG. 5. Self-energy diagrams giving the third-order corrections for the optical magnon gap. The vertices are denoted by their 
superscripts. The diagrams 61, 63, c, dr, and dg represent groups of diagrams which can be obtained by the following vertex 
substitutions: 61) (4,8,7) ^ (6,6,8); 63) (5,2,7) ^ (6,3,8); c) (5,4,+) ^ (6,4,-), (6,8,+) (two diagrams), (5,7,-) (two diagrams); 
d^) (4,8,7) ^ (7,8,4), (7,4,8); ds) (5,6,4) ^ (6,6,8), (5,5,7)(see also the notations in Fig. |). 

The analytic expressions related to the diagrams in Fig. ^ can easily be obtained by means of standard diagrammatic 
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rules (see, e.g., Ref. ^l]). For example, the diagrams oi, di, and ei in Fig. || give the following contributions for the 
dispersion of optical magnons: 



(2S')3 N 



P 



r (a 

[^p 



) , ,,(/3) 



(31) 



(25')3 7V3 



E 

1-5 



T^(7) y(9) y{8) 
^12:34 ''34;5/c'^fc5;12 
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{a} 



- UJ. 



(a) 



,(/3) 



(/3) 



(a) 



,(/3) 



,(/3) 



(32) 



'^^''^(ei) 



1 



(2S')3 7V2 



(+)t/(5) t.(4) 
*^34;2fc '■^2fe;43 



,(/3) 



+ ^^2 



(33) 



The results for ^3 in systems with different site spins (S'l, S2) are collected in Table 1. It is seen that the third-order 
corrections for A are numerically small (as compared to the second-order correction) even for the extreme quantum 
system (^1,52) = (l,i). 



IV. SUMMARY OF THE RESULTS AND 
DISCUSSION 

The spin-wave results for the spin-stiffness constant 
Ps and the magnon gap A for a number of combinations 



5*2) are summarized in Table 1. The restJis-for A are 
compared with available DMRG estimates.EBE^I We find 
an excellent agreement with the numerical estimatesL^I 
for a number of systems, the largest deviation (about 
0.88%) being connected with the system (1, i). 



TABLE I. Spin-wave results for the spin stiffness constant ps and optical magnon gap A of quantum Heisenberg ferrimagnetic 
chains containing two different spins Si > 5*2 in the elementary cell. (ri,r2) and (52,53) are, respectively, the coefficients of 
the spin-wave series for rs = ps/JaoSiS2, Eq. (||), and 5 = A/2(5'i - S2), Eq. @. S' indicates numerical estimatesE3 for the 
reduced excitation gap 5. ri"^ = p^s*^ / JaoSiS2- 



(5l,52) 


Si 


S2 


S3 


S 


s' 


ri 


r2 


^(0) 

' S 


rs 




0.6756 


0.1095 


-0.0107 


1.7744 


1.76 


-0.2391 


0.0283 


0.7609 


0.7892 


(1,1) 


1.0428 


0.4262 


0.0812 


1.6381 


1.63 


-0.4907 


0.0184 


0.7546 


0.7592 


(l-l) 


0.4013 


0.0251 


-0.0047 


1.4217 


1.42 


-0.0959 


0.0160 


0.9041 


0.9201 


(2,1) 


0.6756 


0.1279 


-0.0103 


1.3685 


1.37 


-0.2391 


0.0391 


0.8804 


0.8902 


(2,1) 


0.2861 


0.0095 


-0.0018 


1.2938 


1.29 


-0.0517 


0.0083 


0.9483 


0.9566 



Since in the extreme quantum cases the constructed 
perturbation expansions basically rely on the suggested 
smallness of the quasiparticle interaction V, it is instruc- 
tive to check the self-consistency of the theory by writing 
the series for ps and A in terms of the formal parame- 
ter A = 1, Eq. (pT|). Using the results from Table 1, the 
series for the system (S'l, 6*2) = (1, ^) read 



Ps 



J aoS'iS'2 



0.7609A° -I- 0.0283A2 -I- O (A^) 



(34) 



2J(5i - S2) 



1.6756A° + 0.1095A2 - 0.0107A3 + O (A^) 

(35) 



The above expansions explicitly demonstrate that the 



quasiparticle interaction V in Eq. ( |11| ) introduces nu- 
merically small corrections to the zeroth-order princi- 
ple approximation based on the quadratic quasiparti- 
cle Hamihonian Ho = C(A°), Eq. (|l|). Of course, the 
smallness of corrections by itself does not ensure a good 
quality of the spin- wave expansion. The main weakness 
of the spin-wave theory is connected with the assump- 
tion that the long-range order is well established: it in- 
cludes only transverse spin fluctuations, whereas the lon- 
gitudinal spin fluctuations are completely neglected. In 
this respect, a typical example arises when spin-wave ex- 
pansions are used for magnetic systems near the order- 
disorder critical point where the transverse and longitu- 
dinal fluctuations should be treated on equal ground. In 
spite of the fact that the corrections are small, thCpSpin- 
wave series give unsatisfactory quantitative results. O On 
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the other hand, as the distance from the critical point 
increases the spin-wave description becomes more and 
more reliable. The indicated discrepancy for the acous- 
tical branch in the extreme quantum system (1, i) (al- 
though the second-order corrections to the principle ap- 
proximation are numerically small) probably reflects the 
discussed weakness of the SWT. As a matter of fact, tak- 
ing the parameter Mq as a measure of the distance from 
the disordered phase, such a behavior of the spin-wave 
series in the (Si,S2) family of ID quantum Heisenberg 
ferrimagnets can also be indicated for the magnon gap 
A (see Table 1) and the parameters Eq and mi (see Ref. 
p3| ): the largest deviations from the DMRG results ap- 
pear in the cases with minimal magnetization density, 
Mq = 1/Aao. However, even in the extreme quantum sys- 
tem (1, |) the discrepancies are small and the spin- wave 



expansion produces precise quantitative results. Con- 
cerning the dispersion of acoustical branch w^"^ (and 
the related spin-stiffness constant ps), more numerical 
results are needed to make a statement about the accu- 
racy of spin-wave description in this case. We believe, 
however, that at least for the systems with AIq > l/Aag 
the reported results for ps closely approximate the true 
spin-stiffness constants at zero temperature. 
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APPENDIX A: DYSON-MALEEV VERTICES FOR THE FERRIMAGNETIC MODEL 



Using the symmetric form adopted in Ref. the normal-ordered Dyson-Maleev quasiparticle interaction Vdm, 

J 



smg 
Eq. (y), reads 



Vdm = ^ <^?2 [^i2;34"i"2"3a4 + iV^^^ziC^iha^ai -f 2vl^]^^a\a\0la4, 

1-4 

+ W^%Ao^AP2 + 2Vl%Ao^^l32l3i + 2V^%Ap14Pi 

+ V^^^^Ao^IpIPa + Vi^^piha^ai + V^^^Al^ll^^Pi] ■ (Al) 
The explicit form of the symmetric vertex functions for the fcrrimagnctic Heisenberg model (^ is 

^12-34 = "l"2U3U4"^/2^34 , i=l,2,...,9, (A2) 



where 



^12?34 = + 71-32^12:3 + 71-42:124 -I- 72-32:22:3 + 72-42:22:4 

- o-^^^(7i-3-42:i2;3a;4 + 72-3-42:2230:4) - Cr"^/^(7i2i -I- 7222) , 

*^12;34 = - 71-3212223 - 7i_42i2224 - 72-323 " 72-424 

+ Cr^/^(7i„3_42i222324 + 72-3-42324) "f d"^/^ (71 2i22 + 72) , 

t7(3) 

'^12;34 = - 71-321 - 71-4212324 - 72-322 - 72-4222324 

+ 0-^^^(71-3-42124 + 72-3-42224) + 0-^^/^(712123 -I- 722223) , 
^/2';34 = + 71-321222324 + 71-42122 + 72-32324 + 72-4 

^^^(71-3-4212223 + 72-3-423) - 0-^^/^(71212224 -I- 7224) , 



— a 



^12% = ^ 71-3222324 - 71-422 - 72-3212324 - 72-421 

+ 0-^/^(71-3-42223 -I- 72-3-42123) -f 0-"^/^(7i2224 -I- 722124) 

'^12;34 = - 71-324 - 71-423 - 72-3212224 - 72-4212223 

+ 0-^/^(71-3-4 + 72-3-42122) + 0-^^/^(712324 -I- 7221222324) . 

^12;34 = + 71-32124 + 71-42123 + 72-32224 + 72-42223, 

- 0-^/^(71-3-421 -I- 72-3-422) - 0-^^/^(71212324 -I- 72222324) . 

^1*2-34 = + 71-32223 + 71-42224 + 72-32123 -I- 72-42124 
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- <J^^^{'ji-3-4X2X3X4, + 72-3-4a;ia;32;4) - cr ^/^(7ia;2 + 72^1) , 

~ (9) 

'^12;34 = + 71-32^22:4 + 71-42:22:3 + 72-32:12:4 + 72-42;i2;3 

- f7^/^(7i_3_42;2 + 72-3-42:1) - cr"^/^(7i 2:22:32:4 + 722:12:3X4). 



* Permanent address: Institute of Solid State Physics, Bul- 
garian Academy of Sciences, Tzarigradsko chaussee-72, 
1784 Sofia, Bulgaria. 
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